Modelling the atomic structure of very high-density amorphous ice 
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The structure of very high-density amorphous (VHDA) ice has been modelled by positionally 
disordering three crystalline phases, namely ice IV, VI and XII. These phases were chosen because 
only they are stable or metastable in the region of the ice phase diagram where VHDA ice is 
formed, and their densities are comparable to that of VHDA ice. An excellent fit to the medium 
range of the experimentally observed pair-correlation function g{r) of VHDA ice was obtained by 
introducing disorder into the positions of the H2O molecules, as well as small amounts of molecular 
rotational disorder, disorder in the 0-H bond lengths and disorder in the H-O-H bond angles. The 
low-A; behaviour of the experimental structure factor, S{k), is also very well reproduced by this 
disordered-crystal model. The fraction of each phase present in the best-fit disordered model is very 
close to that observed in the probable crystallization products of VHDA ice. In particular, only 
negligible amounts of ice IV are predicted, in accordance with experimental observation. 

PACS numbers: 61.43.-j 
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The study of water has attracted much attention 
This is not only because of its importance in nature, but 
also because it has a number of unusual properties. In 
particular, it has a very rich phase diagram. In addi- 
tion to liquid water, thirteen crystalline phases of ice are 
known to exist ^s well as at least three amorphous, or 
glassy, phases |3- Mishima et al. ^ were the first to 
synthesize high-density amorphous (HDA) ice, and the 
relationships between this phase, low-density amorphous 
(LDA) ice, and Hquid water have been well studied 0. 
Recently, another distinct phase of amorphous ice has 
been made, by heating HDA ice isobarically at 1 — 2 GPa 
from 77K to ~ 165K. The new phase of ice can be recov- 
ered at 77K and 1 bar |^ . This new phase is even higher 
in density than HDA ice, and has become known as very 
high-density amorphous (VHDA) ice j^. VHDA ice, and 
its relationship to HDA and LDA ice, is of considerable 
interest in studies of polyamorphism [6j (in which a sub- 
stance exhibits more than one amorphous phase), and its 
structure is also relevant to the structure of liquid water. 

Water molecules form strong hydrogen bonds to neigh- 
bouring water molecules, resulting in a roughly tetra- 
hedral molecular arrangement. The hydrogen-bond- 
induced flexibility of this structure is the main reason 
for the many different phases of ice. HDA and VHDA 
ice exhibit closer packing of the water molecules. Each 
water molecule is still bonded to four nearest neighbours 
but, between the first and second neighbour shells, there 
is approximately one extra molecule in HDA ice, and 
two extra molecules in VHDA ice. This increase in co- 
ordination number is observed experimentally 0) IM 
and attributed to occupation of the "lynch pin" location 
by interstitial water molecules which are not hydrogen 
bonded to the origin molecule. Recent MD simulations 
support this interpretation , and suggest that VHDA 



ice is more stable than HDA ice, and that it is the amor- 
phous analogue of the postulated high-density liquid wa- 
ter structure [1^ . 

The medium-range structure of amorphous materials, 
that is how atomic positions are correlated over several 
interatomic distances, is of much interest VHDA 
ice is exceptional among amorphous materials in exhibit- 
ing very pronounced medium-range order (MRO). The 
orientationally averaged real-space pair-correlation func- 
tion, g{r), for r > SA, consists of damped real-space os- 
cillations extending to 20A, with only one significant 
period, D 13] . The orientationally averaged structure 
factor, S{k), measured by neutron and X-ray diffraction 
[l^. exhibits an intense, sharp peak at fco = 2.29(2)A~^ 
in the neutron data, and 2.30(2)A~^ in the X-ray data, 
this peak being the dominant feature (particularly for 
the neutron data). Since g{r) and S(k) are related by a 
Fourier transform, the period of the oscillation in g{r) is 
related to the position of the peak in S{k) hy D — 27r//co- 

The structure of amorphous materials has recently 
been modelled by adding Gaussian positional disorder 
into the atomic structure of crystalline phases of these 
materials [3 ^| • These disordered crystals gave good 
representations of the MRO of some amorphous mate- 
rials, but not necessarily the short-range order (SRO). 
In particular, the value of the wavevector transfer of the 
so-called first sharp diffraction peak (FSDP, a character- 
istic peak in the S{k) of many amorphous materials at 
low wavevector, k ~ 1 — 2A~^) was successfully predicted 
in amorphous silicon and vitreous silica. The existence of 
the FSDP is often regarded as a signature of MRO • 

For a model crystalline lattice, S{k) consists of a se- 
ries of ^-functions caused by scattering from the various 
Bragg planes. The effect of Gaussian positional disorder 
is to multiply these peaks by the Debye- Waller factor. 
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e~ , where a is the half width of the Gaussian distri- 
bution, and to produce a low-amplitude diffuse scattering 
[l^ . It is clear that the high-fc peaks will be destroyed 
first by introducing disorder, and that for large enough 
disorder (high enough a), only one peak (that at the low- 
est value of k) will contribute substantially to S{k). 

Similarly, for a model crystalline lattice, g{r) consists 
of (5-functions at the various interatomic distances. The 
effect of Gaussian positional disorder is to convolve these 
peaks with the Fourier transform of the Debye- Waller 
factor, which is also a Gaussian For large enough 

disorder (high enough a) , these Gaussian peaks broaden 
and combine to form damped oscillations in g{r) with 
a single period [13. This functional behaviour, of 

course, results from the Fourier transformation of S{k) 
having only one significant peak. 

The experimentally measured g{r) for VHDA ice [l3| . 
showing pronounced MRO, characterized by nearly pe- 
riodic oscillations extending to large r, accompanied by 
a dominant sharp low-fc peak in S{k), are very reminis- 
cent of the simulated data for positionally disordered lat- 
tices ^1 • It was natural to ask whether the MRO in 
VHDA ice could be modelled by disordering crystalline 
ice. The structure of an amorphous phase is most likely 
to resemble the crystalline phase with a similar den- 
sity that is stable at similar temperatures and pressures. 
VHDA ice is formed at ^ 165K and 1-2 GPa, with den- 
sity p = 1.25 g/cm"^ (all densities are for H2O) when 
recovered to 1 bar at 77K 0]. According to the phase 
diagram, ice VI is the only crystalline phase stable in 
this region 0. Four other crystalhne phases (Ic, IV, IX 
and XII) are metastable |3| and hence their occurrence 
around these temperatures and pressures cannot in prin- 
ciple be excluded. However, ice Ic and IX can be elimi- 
nated as candidates because their densities (respectively 
0.93 g/cm3 at 77K and 1 bar 0, and 1.17 g/cm^ at llOK 
and 1 bar are so much lower than that of VHDA 
ice, and it is very unlikely that a positionally disordered 
phase could have a larger density than the corresponding 
crystalline phase. If there were significant density varia- 
tions in VHDA ice, then small-angle scattering would be 
experimentally observed, which it is not. It was there- 
fore decided to disorder the structures of ice IV (p = 1.29 
g/cm3 at 260K and 0.5GPa [13), VI (p = 1.37 g/cm^ at 
225K and l.lGPa ^) and XII [p = 1.29 g/cm^ at 260K 
and O.SGPa 0|), in order to try to simulate VHDA ice. 

Crystal models (with ~ 3 x 10'' atoms) of each phase 
were created. The crystalline atomic positions for ice 
IV 0, VI 1131 and XII ^ were taken from neutron- 
diffraction data. The hydrogen atoms were positioned 
according to the "ice rules" |2l|: along each 0-0 bond 
there is only one hydrogen atom, which is covalently 
bonded (and thus closer) to one oxygen atom, and hy- 
drogen bonded to (and thus further away from) the other 
oxygen atom. For these phases, this means that only 
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FIG. 1: The computed pair-correlation function, g{r), for po- 
sitionally disordered ice IV, VI and XII lattices, compared to 
experimental data for VHDA ice. The parameters used to 
generate the simulated curves are those of the best-fit three- 
phase models, as given in Table |I] The dashed vertical fines 
show the peak positions in g{r) for VHDA ice. It is clear that 
the peak positions and periods of oscillation are different for 
each phase. The g{r)'s have been offset vertically for clarity. 



half of the sites available for hydrogen atoms are occu- 
pied; this fractional occupancy is observed experimen- 
tally 0,^3- 

An algorithm due to Buch et al. |22| was 
used to position the hydrogen atoms on their crystalline 
sites. The models were rescaled slightly so that their 
densities matched that of VHDA ice. 

Disorder was then introduced into the models as fol- 
lows. Positional disorder was introduced by displacing 
the H2O molecules rigidly in the x-, g- and z-directions. 
Then the H2O molecules were rotated around the x-, g- 
and z-axes passing through the oxygen atom positions, 
such that the oxygen atoms were kept stationary. The 
hydrogen atoms were then moved along the 0-H bond 
direction, to disorder the 0-H bond lengths. Finally, the 
hydrogen atoms were symmetrically moved in the H-0- 
H plane to disorder the H-O-H bond angles. In each 
case, the change in the relevant quantity (position or an- 
gle) was randomly chosen from a Gaussian distribution 
centred on zero, with half widths as given in Table 

Calculated g{rys for these models were fitted to 
neutron-diffraction data of deuterated VHDA ice, 
the structure of which differs negligibly from the struc- 
ture of the hydrogenated material 23] . The neutron data 
were used, as they were of a higher quality than the X-ray 
data. Since particularly good agreement with the data 
at low r was not expected, the data for r < 7 A were ex- 
cluded from the fit. The partial pair-correlation functions 
900 gooif) and (r) of our models were combined 
to give the full neutron g{r) as g{r) = 0.088g 00 {'>')+ 
0A18goDir)+ OAMgooir), as in Ref. 0. 

The calculated (?(r)'s for the positionally disordered 
models of each separate phase do not exactly fit the g{r) 
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Half width 


Range studied 


ice IV 


ice 


VI 


ice 


XII 








3 phases 


2 phases 


3 phases 


2 phases 


3 phases 


Positional disorder 


^pos 


U.Uo — U.ouA 


0.22A 


0.16A 


0.16A 


0.19A 


0.18A 


Rotational disorder 




0.01 - 0.05 rad 


0.04 rad 


0.04 rad 


0.04 rad 


0.03 rad 


0.01 rad 


O-H bond length disorder 


<JOH 


0.01 - 0.04A 


0.03A 


0.04A 


0.04A 


0.02A 


O.OlA 


H-O-H bond angle disorder 


choh 


0.01 - 0.04 rad 


0.01 rad 


0.01 rad 


0.01 rad 


0.03 rad 


0.02 rad 


Volume fraction, / 




0-1 


0.035 


0.544 


0.549 


0.456 


0.416 



TABLE I: Disorder parameters for the crystalline phases, ice IV, VI and XII. The ranges of values of the disorder parameters 
and the fractions of each phase, and the values in the best-fit two- and three-phase models, are also given. In each case, the 
change in the position or angle was randomly taken from a Gaussian distribution, centred on zero, with half widths as shown. 



of VHDA ice (see Fig.^. 

Hence, the structure of VHDA ice was assumed to com- 
prise an inhomogeneous mixture of regions of position- 
ally disordered ice IV, VI and XII, of indeterminate size, 
but with total volume fractions of fiv , fvi a-nd fxii = 
1 ^ ifiv + fvi), respectively. The g{r)'s of the three 
phases were combined to give a total pair-correlation 
function: g(r) = fivgiv{r)+ fvigvi(r)+ fxiigxii{r), 
where < fiv, fvi, fxii < 1, and e.g. giv{r) is the full 
g{r) of the disordered ice IV lattice. For each combi- 
nation of the g{r)'s of the three phases, a least-squares 
minimization was performed with respect to fiv and fvi- 

The minimization analysis showed that dpos and the 
volume fractions were the most significant parameters in 
determining the quality of the fits in the medium range 
(r > 7 A). The other parameters, ag, aoH and cthoh, 
were much less important, and only significantly affected 
g{r) at low r. The best-fit model (with parameters as 
given m Table P) gives a total g(r) with medium-range 
oscillations in very good agreement with those found ex- 
perimentally for VHDA ice (Fig. |2J), and very signifi- 
cantly better than for the disordered individual phases 
(see Fig. ^ . The short-range fit is less good; in partic- 
ular, the heights of the peaks at lA and ~ 3.3A, are 
not reproduced so well (see inset to Fig. [21). The height 
of the latter peak is known to be sample and pressure 
dependent, so one would not necessarily expect a good 
fit. The best-fit model has a least-squares deviation from 
the experimental data of 0.0777. It contains 3.5% ice IV, 
an amount which is not resolvable by diffraction studies. 
The fitting procedure was repeated for disordered ice VI 
and XII only, and the best-fit model (Fig. ^ has a devi- 
ation of 0.0779, which is only 0.25% more than the devi- 
ation for the best three-phase model. The corresponding 
g{r)'s are indistinguishable (Fig. [^J. We conclude that 
the amount of ice IV needed to model VHDA ice is negli- 
gible, and that disordering only ice VI and XII produces 
just as good a description of the MRO of VHDA ice. 

It is also instructive to compare the structure fac- 
tor S{k) of the best-fit model (calculated by Fourier- 
inverting g{r)) with the experimental data (Fig.OJ. The 
height, width and position of the first peak are fitted 
very well, which is expected since this peak governs the 
medium-range real-space behaviour. The fit at higher k 
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Best-fit Iwo-ptiase model 

— Best-fit tliree-phase model 

— Experimental VHDA ice 




FIG. 2: The pair-correlation functions g{r) of the best-fit 
models consisting of disordered lattices of ice IV, VI and XII 
(dashed line) and ice VI and XII only (dotted line), com- 
pared with the experimentally measured neutron-diffraction 
g{r) of VHDA ice [l^ (solid line). The disorder parameters 
of the models are given in Table |l] The best-fit two- and 
three-phase models are indistinguishable for medium-range 
distances. The g{r) for VHDA ice was calculated by Fourier 
inversion of the experimentally measured structure factor [1^ 
without smoothing, resulting in considerable noise, and spu- 
rious peaks for r < O.SA. 



(which corresponds to small r) is less good, although the 
peaks are reasonably well reproduced. 

When Klotz et al. recrystallized high-pressure HDA 
(probably VHDA) ice, the fractions of the phases present 
were 56% ice VI and 44% ice XII These are very 

close to the simulated values of 54% ice VI and 46% ice 
XII for the best-fit two-phase model. This implies that 
the structure of VHDA ice is intimately related to the 
structures of ice VI and XII. It is also interesting to note 
that the value of apos in our models (Table P) is compa- 
rable to the thermal mean-square atomic displacement, 
which for these phases is ^ 0.2 A 19., ..20j . 

Attempts were made to improve the fit at short range, 
by fitting separately to the low-r region of g{r), and to 
the medium- and high-fc region in S'(fc), whilst varying 
the short-range disorder parameters ag , uqh , cfhoh and 
volume fractions, but none of these attempts improved 
the fit of the model. Some of this discrepancy may be 
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FIG. 3: The structure factor S{k) of the best-fit model con- 
sisting of disordered lattices of ice IV, VI and XII, compared 
with the ex peri mentally measured neutron-diffraction S{k) of 
VHDA ice [ij]. The parameters of the model are given in 
Table El 



attributed to effects at the boundaries between regions of 
disordered ice phases proposed in the VHDA ice struc- 
ture. Our models take no account of these boundaries, 
yet if the ice structure consists of positionally disordered 
domains of different crystal phases, then the interfacial 
regions will not necessarily resemble any of the phases. It 
is likely that any such distortion will be fairly local, ex- 
tending over maybe two nearest- neighbour bond lengths, 
confining these surface effects to low r. 

The fact that the disordered crystal model is able to 
reproduce the medium-range structure, but not partic- 
ularly that of the first- and second-peak region in g{r), 
for VHDA ice may be an important indication of struc- 
tural differences between the amorphous and crystalline 
phases. It is well known that the 4-8 A range in amor- 
phous ice is very sensitive to changes in pressure, and X- 
ray diffraction studies of VHDA ice have shown that this 
region contains two peaks which are most likely related to 
the start of the formation of non-bonded interpenetrating 
lattices (similar to those present in ice VI and XII). A dis- 
tortion of the local tetrahedral unit in VHDA ice enables 
an increased packing of interstitial molecules, suggesting 
that much of the pressure-induced frustration in the ma- 
terial occurs largely in the ordering of the molecules just 
beyond the first shell. This frustration arises from com- 
petition between the optimum local packing of molecules 
and the preferred longer-range topology of the network. 
Such high-density structures are very different from that 
characterized by the well-defined first two shells found in 
LDA ice, which closely resembles that of ice Ih. 

In conclusion, the medium-range structure of very 
high-density amorphous (VHDA) ice has been modelled 
by adding positional disorder to three crystalline phases 
of ice, namely ice IV, ice VI and ice XII. These phases 
were chosen because they are stable or metastable in the 



region of the phase diagram where VHDA ice is formed, 
and their densities are similar to that of VHDA ice. For 
medium-range distances (r > 7A), the fit to the experi- 
mentally observed neutron-diffraction g(r) is very good. 
The fractions of each phase present in the model which 
gives the best fit to the experimental data are very close 
to those found when a sample, which was likely to have 
been VHDA ice, was recrystallized in particular, 

the amount of ice IV present is negligible (none was ob- 
served experimentally). We conclude that the structure 
of VHDA ice at large enough distances from an arbitrary 
origin atom can be well modelled as an inhomogeneous 
mixture of positionally disordered ice VI and ice XII. 

JKC would like to thank the Engineering and Physical 
Sciences Research Council for financial support. 
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